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We analytically study rogue-wave (RW) solutions and rational solitons of an integrable fifth-order 
nonlinear Schrodinger (FONLS) equation with three free parameters. It includes, as particular cases, 
the usual NLS, Hirota, and Lakshmanan-Porsezian-Daniel (LPD) equations. We present continuous- 
wave (CW) solutions and conditions for their modulation instability in the framework of this model. 
Applying the Darboux transformation to the CW input, novel first- and second-order RW solutions 
of the FONLS equation are analytically found. In particular, trajectories of motion of peaks and 
depressions of profiles of the first- and second-order RWs are produced by means of analytical and 
numerical methods. The solutions also include newly found rational and W-shaped one- and two- 
soliton modes. The results predict the corresponding dynamical phenomena in extended models of 
nonlinear fiber optics and other physically relevant integrable systems. 


Solitons are usually generated by the balance 
of linear dispersion and nonlinear self-focusing. 
Rogue waves (RWs), as a special type of solitary 
waves, are driven by a similar mechanism act¬ 
ing on top of a modulationally unstable flat back¬ 
ground (alias continuous wave, CW). In particu¬ 
lar, the celebrated nonlinear Schrodinger (NLS) 
equation with the self-focusing sign of the non¬ 
linearity gives rise to both bright solitons and 
RWs, which are used to model nonlinear phe¬ 
nomena in diverse flelds, such as nonlinear op¬ 
tics, deep ocean, plasmas, Bose-Einstein conden¬ 
sates, and biophysics. In nonlinear flber optics, 
the correct description of tall narrow pulses pro¬ 
duced by the RWs makes necessary to include 
higher-order effects, such as the stimulated Ra¬ 
man scattering, nonlinear group-velocity disper¬ 
sion (alias the shock term), and third-order dis¬ 
persion. In the general case, the addition of such 
terms to the NLS equation breaks its integra- 
bility, making it impossible to And solitons and 
RWs in an analytical form. However, under spe¬ 
cial conditions the respective extended NLS equa¬ 
tion may keep the integrability. In particular, 
quite a general model of the latter type, namely, 
the fifth-order NLS (FONLS) equation, which in¬ 
cludes three free parameters, was recently found. 
Its particular forms amount to other previously 
known integrable equations, such as the usual 
NLS, Hirota, and Lakshmanan-Porsezian-Daniel 
(LPD) equations. In this paper, our objective is 
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to use the integrability of the FONLS equation 
for constructing its exact RW solutions. Because 
RWs are driven by the modulational instability of 
flat CW states, we first present the CW solutions 
of this model and conditions for their modulation 
instability. Then we proceed to the core part of 
the analysis, which reveals novel exact first- and 
second-order RW solutions, by means of the Dar¬ 
boux transformation applied to the CW input. 
The solutions are obtained in the form of explicit 
rational functions of the spatial coordinate and 
time. In particular, interesting results are pro¬ 
duced by the analysis of trajectories of the motion 
of peaks and depressions of the RW profiles. Fur¬ 
thermore, the obtained rational solutions include 
new W-shaped one- and two-soliton modes of the 
FONLS equation. The analysis developed in this 
work and results produced by it may suggest new 
possibilities to analyze other sophisticated non¬ 
linear models of practical significance. 

I. INTRODUCTION 

Rogue waves (RWs, also known as monster waves, 
killer waves, giant waves, extreme waves, etc.) 0 , 0 , as 
a kind of localized modes, are generated, in the paradig¬ 
matic form 0 , by the self-focusing nonlinear Schrodinger 
(NLS) equation for wave amplitude {x,t), 

= 0 . ( 1 ) 

as a limit case of the Ma-breather Q or Akhmediev- 
breather Q solutions on top of a flat continuous-wave 
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(CW) state, which is subject to the modulational insta¬ 
bility (MI) [5[. This equation is written in terms of the 
guided-wave notation, with x and t being thy^ropagation 
distance and reduced time, respectively i0. The RWs 
are sometimes called rogons^ if they reappear unaffected 
in the size and shape after their interactions [H. 

The RW phenomenon was initially found in the deep 
ocean [l|, [H . Recently, it has also drawn growing interest 
in other fields, such as Bose-Einstein condensates (BEC), 
nonlinear optics, plasma physics, and even dynamics of 
financial markets (Sin. Generally speaking, the RWs 
are spatially and temporally localized excitations which 
seem “to app ear from nowhere and disappear without 
a trace” [l 6 |, that makes it difficult to systematically 
observe them in the ocean, where they had been origi¬ 
nally spotted. The origin of the RWs is still a subject 
of discussions 0 , stimulating the development of var¬ 
ious theoretical approaches and laboratory experiments 
[m, [HI. In particular, it is well established that the MI 
of the CW background, also known as the Benjamin-Eeir 
instability | 2 Q||. plays a pivotal role in the generation of 
RWs 

In addition to the basic RW solutions of the NLS 
equation, RWs have been found in related integrable 
equations including higher-order dispersion and nonlin¬ 
earity terms, such as Hirota [HI, Sasa-Satsuma [H, 
Ger dj ikov-Ivanov [iilsil , Lakshmanan-Porsezian- 
Daniel (LPD) [H|, quintic NLS [30|, modified NLS 
im, and Davey-Stewart son [HI equations, the deriva¬ 
tive NLS equations of the Chen-Lee-Liu and Kaup-Newell 
types [HI, coupled NLS equations 35, [HI, three-wave 
equations [s^, and various modifications of these mod¬ 
els with variable coefficient 38-4o|. Einding RWs and 
rational solitons of other nonlinear wave equations in the 
field of nonlinear science remains a relevant subject. 

Recently, Chowdury et al 0 reported an inte¬ 
grable three-parameter fifth-order NLS (EONLS) equa¬ 
tion, written as 


+ S{2p) - iaH{2p) -h 7^('0) “ = 0, (2) 

where ^ is a complex field, subscripts denote 

partial derivatives with respect to the corresponding vari¬ 
ables, while (a, 7 , and 6 are all real parameters. Eurther, 
S{ 2 p) denotes the NLS part 

Siip) = ^V’tt + IV’PV’, 

iL('0) stands for the Hirota part, 

H{7p) = 

Pi'ip) is the LPD part. 


and Q(' 0 ) is the quintic part. 


QW = ^ttttt + + 10 1 ^ 7 )^ 

where t/i* stands for the complex conjugate of This 
equation 0 is a member of the NLS hierarchy intro¬ 
duced by Kano 0 and contains many special cases. 
When (a = 7 = ^ = 0, Eq. m reduces to the funda¬ 
mental NLS equation {p. As 7 = (5 = 0 , it becomes 
the Hirota equation [^. When a = (5 = 0, it is the 
LPD equation [i^. All these particular nonlinear wave 
equations are integrable, and the respective soliton solu¬ 
tions have been obtained using different methods. Equa¬ 
tion ([2]) with Of = 7 = 0 is also integrable, and its soliton 
and breather solutions are available too [30|, 1^ . 

To the best of our knowledge, rational RW solutions of 
Eq. ([ 2 ]) with arbitrary parameters (a, 7 , and (5) have not 
been reported before. In this work, we aim to produce 
such localized solutions and study their dynamical prop¬ 
erties by means of the generalized Darboux transform. 

The rest of the paper is organized as follows. In Sec. 
H, the Lax pair and Darboux transform with nonzero 
CW background are introduced for Eq. ([ 2 |). In Sec. HI, 
the MI of the CW states is investigated, as the driving 
mechanism behind the existence of RWs. In Sec. IV, 
the first- and second-order RW solutions for Eq. m are 
explicitly found by means of the Darboux transform ap¬ 
plied to the flat CW. Trajectories of motion of peaks and 
depressions of the derived RWs are explicitly considered, 
to demonstrate the structure and dynamics of these pat¬ 
terns. Rational W-shaped solitons of the EONLS equa¬ 
tion, which differ from the usual bright and dark solitons, 
are also found for special parameter values. The paper is 
concluded by Sec. V. 


II. THE LAX PAIR, DARBOUX TRANSFORM, 
AND CW SOLUTIONS 


Equation ([2]) is integrable in terms of the correspond¬ 
ing Lax pair, which can be obtained as the compatibility 
condition, R^t = Rtx, i-e.,, - Vt + [V, V] = 0 , for the 

following linear spectral problem |4l| : 

Rt = UR, U = XJ + Uo, (3) 

Rx = VR, V = XU + Vo^aL^-fM ^ SN, (4) 

where A is a complex spectral parameter, the dost func¬ 
tion R = R(x, t) is a column vector given by 

/ r \ 


P{ip) = + 6|V’|7 + 


R = R(x, t) = 


s 


(5) 
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with r = r{x,t) and s = s{x,t) being complex functions, 
and matrices J, Uq, Vq, I/, M, N given by 


the Darboux transform ([8j) for obtaining RW solutions of 
Eq. m, we here take the CW solution of Eq. 


j = 
Uo = 


i 0 
0 -i 


Vo = 


where 


in = 


Mo = 


No = 


0 #* 

iip 0 

1 ' V’t 

2 _ -'ipt ijV’P 
L = -4(AV + X^Uo + XVo) + Lo, 
J\d = ‘2XL “t" 

N = -2XM + No, 


V’V’t - i (2|'^| V* + 'ipu) 


mi m2 
ms —mi 

ni n2 
ns -ni 


mi = - i (3|'^|^ + ipiptt - 'ijjt'il’t + 

m2 = — ms 

= 6 iV’pV’r +'>Pttt^ 

m ={tp - ip*)ipttt - ipftl’u + 

n2 = -nl 


(6) 


(7) 


=6i|V’|V* + + 4i|'^trV’* + hiV’V’i 

+ 8 i|^| Vtt + ii’tttf 


/,*2 


Based on the Lax pair of an integrable nonlinear wave 
equation, the Darboux transform is a powerful method 
for generating soliton and breather solutions, using a zero 
solution or a flat CW as the “seed” solution [^. The 
hierarchy of multisoliton states [iHl can be derived start¬ 
ing from zero, while solutions related to the MI are 
obtained starting from the CW. 

The Darboux transform for Eq. m is given by 


'tpn+l{x,t) = i^n{x,t) - 


4isn+i{x,t)r*_^_^{x,t) 
\rn+l{x,t)\^ + \Sn+l{x,t)\'^ 


( 8 ) 

where n = 0,1,2,..., is a seed solution of 

Eq. ([2]), while V^n+i(^, t) is a new solution generated from 
while rn+i(x,t) and Sn-\-i{x,t) are solutions of 
the linear spectral problem based on Eqs. with 

'0(x,t) substituted by We here choose X = i 

in system as this value of the spectral parameter 

corresponds to RW solutions of Eq. ©■ 

The zero seed solution was employed to build multi- 
soliton solutions of Eq. ([2]) in Ref. . In order to use 


Mx,t) = (9) 

as the seed, where a 7 ^ 0 and c are arbitrary real con¬ 
stants, and 

b = ac ( 6 a^ — c^) -h 7 ( 6 a^ — 12 a^c^ + c^) 

^2 (10) 

-\-cS (SOa"^ — 20 a^c^ + ^ 

Lastly, we note that EONLS equation (| 2 ]) is invariant 
under the following scale transformation: 

9 / 

X —y fi x^ t — jjjt^ Ip —y —, 

o, —y fjjd^ 7 —^ l-^^•) ^ ^ 

where ji is an arbitrary real constant. We stress that 
this invariance is different from that for the NLS equa¬ 
tion [ 3 ^. 


III. THE MODULATIONAL INSTABILITY (MI) 
OF CW STATES 


To investigate the MI of the CW solution (|9]), subject 
to constraint m, we add perturbations (p{x,t) with an 
infinitesimal amplitude 5 , 


{p{x,t) = [a-h 50 (x, ( 11 ) 

which is followed by the derivation of the linearized equa¬ 
tion for the perturbation: 

i(px — ^^(pttttt + (7 + ^cS)(ptttt + i( 10 c^ — lOa^^ 

-h4c7 — a)(pttt + (40a^C(5 -h 3ca — 60^7 — lOc^^ 
-l- 8 a ^7 -|- l/2)(ptt + 2a^(5c^ + 7 ) 0 ?t T i(—30a^ 
-\-60a^c^S — 5c^S + 24a^c7 — 4c^7 — 6a^a 
-\-?>c^a -h c)(pt -h a^(60a^C(5 — 20c^S + 12 a ^7 
— 12 c ^7 -h 6ca l){(p + 0*) = 0. 

( 12 ) 

Relevant solutions to Eq. ([ 12 ]) with wavenumber k and 
frequency 00 are sought for as 

0(x, t) = F cos{kx — ujt) + iG sm{kx — oot), (13) 

where F and G are real amplitudes. The substitution of 
ansatz ([13]) into Eq. ([T^ readily leads to the following 
dispersion relation for the perturbations, obtained as the 
condition for the existence of nontrivial solutions for F 
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and G: 

4 {/c — cj [c + a{3c^ — 6 a^ + k‘^) + 47 c( 6 a^ — — k‘^) 

+5(5(12a^c^ — 6 a^ + 2oF‘k‘^ — 2c^uj‘^ — — uj^)] 

= —Cc;^(4a^ — [lOc(5(6a^ — 2 c^ — cj^) 

+ 27 ( 6 a^ — 6 c^ — cj^) + Gac + l] ^ • 

(14) 

It follows! rom Eq. m that we know the wavenumber k 
is complex for frequencies interval 


q/ = (5 = 7 = 0 , we also find the known first-order RW 
solution of the self-focusing NLS equation, 


'01,NLS 


4(1 -f 2ix) 
4(t^ -h x‘^) -h 1 


(ii) for ^ = 7 = 0 , we obtain the first-order RW solution 
of the Hirota equation, 


V^ljHirota 


4(1 + 2ix) 

4(t -h 6ax)^ -h 4x^ -h 1 


|w| < 2\a\, 

hence the CW states are subject to the MI, and RW 
patterns may be expected as solutions of Eq. (ED in this 
case. 


IV. ROGUE WAVES, RATIONAL SOLITONS, 
AND THE DYNAMICAL ANALYSIS 


(hi) for a = h = 0, the first-order RW solution of the 
LPD equation is obtained: 


'01,LPD 


4[1 + 2i(127 + l)a;] ] 

4^2 + 4(127 + l)2a;2 + 1 J 


In the following we will analyze the solution (p! 8 ]) for 
two different cases 7 7 ^ — 1/12 and 7 = — 1 / 12 , which 
display the distinguish profiles of solution (p! 8 |) . 


A. First-order rational solutions and trajectories of 
extreme points 


1. The first-order RW solution with 7 7 ^ — 1/12 


We choose a = 1 and c = 0 in the seed solution given 
by Eqs. © and m, i-e., 

(15) 

Solving the linear spectral problem with X = i 

and the seed solution m yields eigenfunctions 

ri(x, t) = —co[2(t -T vx T- iBx) — 

(16) 

si(x, t) = co[2(t + px -h iBx) -h 

where cq is an arbitrary nonzero constant, and p, B are 
given by 


5 = 127 - 1 - 1 , v = 6{a^56). (17) 

Then, substituting solution (USD and eigenfunctions 
m into the Darboux transform (|8]), we obtain the first- 
order rational solution of Eq. ([2]) as 


With 7 7 ^ —1/12, the rational solution (p! 8 |) of Eq. (j2j) 
is just its the first-order RW solution. To analyze dy¬ 
namical properties of the RW solution, we consider its 
intensity. 




8 — 32(t + px)^ + 325^x^ 
[1 -h 4(t + vx)‘^ + 45^x^]^ 


( 20 ) 


It is easy to see that expression (|20D has three critical 
points. 


(a;i,ti) = ( 0 , 0 ), (3^2,3, ^ 2 ,3) 



which are three real roots of system {9 (^\ipi{x, t )7 /dt = 
0 , d /dx = 0 } (see Appendix A). Eurther, 

the critical point (xifii) is a maximum (peak), as 


52(|^l(x,t)|2) 

dt^ 


-192 < 0, 


ipi{x,t) 



4(1 + 2iBx) 
Di{x,t) 


exp 


-iB + l)x 


(18) 


where we define 


D\ (x, t) — 4(t “h vx)‘^ -|- AB^x^ 1 5 (1^) 

Erom the first-order rational solution (jlSp of Eq. 0 , 
the corresponding first-order RW solutions can be ob¬ 
tained also for the Hirota, LPD, and quintic NLS equa¬ 
tions, by choosing the respective values of a, 7 , and 
(5, as specified in Introduction. Eor example, for (i) 


and 

/ d'^ (|V>i(a;,t)7 Y d'^ (|^^l(x,t)7 d'^ (|^^l(a;,t)7 

^ dxdt J dB dx‘^ 

= - 12288 ( 127 + 1)2 <0 ( 21 ) 

at point (xi,ti), which makes the meaning of condi¬ 
tion 7 7 ^ — 1/12 obvious. Similarly, two critical points 
(^ 2 , 352 , 3 ) are minima (depressions), as 

(|V’i(x,t)|2) _ 

df2 
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FIG. 1. (Color online). The first-order rogue wave (ITS)) for 
a = 7 = (5 = 1: (a) The motion of the peak and depression 
centers, (b) The evolution of the peak absolute value, (c) The 
evolution of the absolute value of the field at the depression 
points. 


2 - 


(a) 

M A 


t -10^ 


■ 0.4 
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FIG. 3. (Color online), (a) The density evolution and (b) 
2D density contour plots for the first-order RW solution (ITSl) 
dX a — ^ — 5 — 1. The corresponding motion of peaks and 
depressions is displayed in Fig. [1] 




FIG. 4. (Color online), (a) The density evolution and (b) 
2D density contour plots for the first-order RW solution (p^ 
at a = —5, 7 = 3 and (5=1. The corresponding motion of 
peaks and depressions is displayed in Fig. [2l 


The evolution of the first-order RW solution can 
be characterized by trajectories of peaks and depres- 


47j . The trajectory of the peak’s center is given 


by the dependence of its temporal coordinate, T/j,, on x. 
It follows from Eq. ([T 8 |) that 


Th = -vx, (23) 

and the temporal coordinates, Tc±, of centers of the two 
depressions are given by 

Tc± = -vx ± 173(1 + 452x2). (24) 


FIG. 2 . (Color online). The same as in Fig. [U but for 
parameters (a = — 5, 7 = 3, and (5 = 1. 

and 

/ ^2 (|V;i(x,()| 2) T g2 (|^^(a;^^)|2) q 2 (|^^ |2) 

y dxdt J dt‘^ dx‘^ 

= -48(127+ 1)^ < 0 (22) 

at points (^ 2 , 3 ,^ 2 , 3 )- 

Note that these critical points, {xj^tj), j = 1,2,3, do 
not depend on parameters (a, 7 , and S. 


Thus, the motion of the peak is totally determined by pa¬ 
rameter V = 6 ((a + 5J) defined in Eq. (fT7|) . while the RW 
width, which is defined as the temporal distance between 
the two depressions, i.e., 

Td= 73(1+452x2), 

is determined by parameter B = I 27 + 1, which is also 
defined by Eq. dm)- Thus, from the full set of three 
constants a, 7 , and 6 in Eq. m, two parameters, a and 
(5, control the motion of the RW peak in solution ([TS]) . 
while the third parameter, 7 , determines the temporal 
distance between the two depressions. 
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FIG. 5. (a) The evolution of the W-shaped profile and (b) 2 D 
density contour plot for rational soliton (|25l) with a = 1 , 7 = 
- 1/12 and (5=1. 


Rational soliton feature W-shaped profiles at all x, 
as shown in Figs. [5] and [6l 

Next, we analyze critical points of the intensity of the 
rational solution (|25]) . It follows from Eq. (|25]) that there 
are three critical points of 

{xi,ii) = {x, -vx), (x 2 . 3 ,( 2 . 3 )= -vx ± , 

(26) 

whose trajectories are straight lines. It also follows from 
Eq. (|25]) that 



fd'^ ill’ll (|V>i-y(a;,()p) ^ 

y dxdt J dt‘^ dx^ 

at three critical points (Ei) [in agreement with Eqs. m 
and ([ 22 ])], hence the points cannot be characterized on 
the basis of this equation. 

On the other hand, the critical point (xi, ti) is a max¬ 
imum (peak), as 


FIG. 6 . (a) The evolution of the W-shaped intensity profile 
and (b) 2D intensity contour plot for rational soliton ()25l) at 
a = —5, 7 = — 1/12 and (5=1. 


To illustrate the structure and dynamical proper¬ 
ties of the first-order RW solution m at different 
parameters, we choose sets ((a = 1,7 = 1 ,^ = !) and 
((a = 5, 7 = 3, (5 = 1) to display traces of peak (|23|) and 
depressions ([M|) in Figs. [T][a) and| 2 ][a). Further, the evo¬ 
lution of the absolute value of the field at the peak and 
depression points is shown in Figs. mb,c) and| 2 ][b,c). 

Lastly, 2 D density contour plots of the RW solution 
(fT8|) are presented in Figs. [3] and ID for different parame¬ 
ters. 


2. The first-order rational soliton with 7 = —1/12 


For the case of 7 = —1/12, we find that rational solu¬ 
tion (fT 8 |) becomes a solitary-wave solution of Eq. m, 


1 pl-y{x,t) 



4 

1 + 4(( + vx)^ 




(25) 


which is a rational soliton-like solution of (EJ with prop- 
erties 


|'0i^(x,t)p ^ 1 as \t-\-vx\^oo^ 
\f;ij{x,t)f ^ 9 as t + vx^O^ 
max(|'0i^(x, t)p) = 9 as t-\-vx = 0^ 


min(|'(/^i^(x, t)p) = 0 as \t + vx\ = V3/2, 


5(2 


-192 < 0 


at this point (see Eigs. [5] and | 6 |). Similarly, two critical 
points (^ 2 , 3 ,^ 2 , 3 ) are minima (depressions), as 


52 _ 

5(2 

at these points (see Eigs. |5]and|6|). 


B. Second-order rational solutions and trajectories 
of extreme points 


To produce second-order RW solution of Eq. (| 2 |), we 
choose the first-order rational solution m found above 
as a new seed. Substituting it into the linear spectral 
problem based on Eqs. ©-© with A = i, we solve it to 
generate new complex functions r 2 (x,t) and 52 (x,t), 

r2ix,t) = [r2i(a;,()-h (r22(a;,()] 

(27) 

S2{x,t) = [S2i(a;,()-|-is22(a;,()] 

where r 2 j{x,t) and S 2 j{x,t) (j = 1 , 2 ) are given by 


r2i{x,t) 


r22{x,t) 


S2l{x,t) 


S22{x,t) 


ci^{x,t) -h C2Bri{x,t) 
4Di{x,t) 

-4c2B^{x, t) -h cir]{x, () 
8Di{x,t) ’ 

4c2i?^(—a;, —() — cir]{—x, —() 
8Di{x, t) 

cii{-x, -() - 1 - C2Bri{-x, -() 
4Di{x,t) ’ 


( 28 ) 
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FIG. 7. (color online). The second-order rogue wave solution 
( [SO]) for a = 2, 7 = 3 and (5 = — 1. (a): The motion of the 
peak and depressions (red and green lines, respectively), (b) 
and (c): The evolution of the absolute values corresponding, 
severally, to the peak and depressions. 


FIG. 8 . (color online). The same as in Fig. [71 but for 
parameters ol — —5, 7 = 3, and (5=1. 


1. The second-order RW solution with 7 7 ^ — 1/12 


^(x, t) = 16Bv{B^ + + 16B{B^ + 3v^)x^t 

-IGB^x^ + ASBvxH‘^ + IGBxt^ 

+4(1205(5 — 3Bv + 4:v)x‘^ 

+(8-205)x + 4(4-75)xt, 
r][x,t) = 16t^{l — t — Avx) ^ 4:S{B‘^v‘^)x‘^t (29) 

+12t — 96v‘^x‘^t^ + 32(30^ + v)xt 
+4(7^ + 120(5)x + 16(5^ - v^)x^ 

+16i;(35^ + v‘^)x^ 

+ 8(752 + 120(5 + 4 ^ 2 - 45 ), 


with expressions 5, 5i(x,t) given by Eqs. ([ITj) and 

Substituting expressions (|28]) and (|29]) into system ([27|) 
and using Darboux transform ([ 8 ]) yields exact second- 
order rational solutions of Eq. m as 


exp Q(H-B)a;^ . 

(30) 

Cumbersome expressions for real functions 52 (x,t), 
G 2 (x,t), and K 2 {x,t) are given in the Appendix B. Par¬ 
ticularly, for cases (i) (a = (5 = 7 = 0 , (ii) (5 = 7 = 0 , 
and (hi) (a = 7 = 0 , solution (HOI yields the second- 
order RW solutions of NFS, Hirota, and LPD equations, 
respectively. 


'tp2{x,t) = 


G2{x,t) +iK2{x,t) 
D2{x,t) 


To analyze dynamical properties of the second-order 
RW solution (pn|i with 7 7 ^ — 1 / 12 , we first consider its 
intensity in the form of 


\'lp2{x,t)\‘^ = 


[D2{x,t) + G2{x,t)]‘^ + K2{x,t) 

Dl(x,t) 


(31) 


Unlike the first-order rational solution ([18]) . the peak 
and depressions of |'(/52(x,t)P cannot be obtained in 
an analytical form for expression m since we need 
to analytically solve two algebraic equations of the 
eleventh order for x and t, which are generated by 
{d{\^ 2 {x,t)\‘^)/dx = 0, d\ ('02(^5))P) /dt = 0}. There¬ 
fore, similar to what was done above for the first-order 
rational solution, we numerically analyze the respective 
trajectories with 7 7 ^ —1/12 (e.g., 7 = 3 ). The results 
are presented in Eigs.[7|and|8l and the corresponding 2D 
contour plots are displayed in Eigs. [QlandfTOl 

A characteristic difference of the second-order RW pat¬ 
tern from its first-order counterpart is that, in a limited 
interval of the evolution variable [namely, \x\ < 0.149 and 
\x\ < 0.152 in Eigs. [3[a) and| 8 ](a), respectively], in addi¬ 
tion to the set of the single peak and two depressions, a 
new peak-depression pair appears. 

To better represent the relevant solutions close to 
X = 0, trajectories of the peaks and depressions, and the 
evolution of the respective absolute values in the interval 
of \x\ < 0.01 are displayed, for different parameters, in 
Eigs. Hi] and nil In the case of (a = 2, 7 = 3 and (5 = — 1, 
there are two peaks and three depressions in the inter¬ 
vals of —0.149 < X < —0.001 and 0.001 < x < 0.149, and 
there are three peaks and four depressions in the interval 
of \x\ < 0.001. In the case of a = —5, 7 = 3 and ^ = 1, 
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FIG. 9. (color online). The second-order rogue wave solution 
([30]) at a = 2, 7 = 3, and (5 = — 1. (a): The density evolution; 
(c) the corresponding 2D contour plot. 



FIG. 10. (color online). The second-order rogue wave so¬ 
lution ( (301) at a = —5, 7 = 3, and (5=1. (a): The density 
evolution; (c) the corresponding 2D contour plot. 


FIG. 12 . (Golor online). The second-order rogue wave solu¬ 
tion (l30)) for a = —5, 7 = 3, and (5=1. (a): The motion of 
the peaks and depressions (red and green lines, respectively); 
(b) The evolution of the absolute values of the peaks and 
depressions (red and green lines, respectively). 

Cumbersome expressions for real functions 
G 2 j{x,t), and K 2 j{x,t) are given in the Appendix C. 

Some profiles of two-soliton solutions (|32|) are displayed 
in Figs. [13] and [14] for (a = 1 , (5 = 1 and a = —5, S = 
which display elastic interactions of two solitary waves. 
It follows from Fig. [14] that the amplitude of one soli¬ 
tary wave grows to a high limit value, and for the other 
solitary wave the amplitude falls to a low limit value at 
\x\ oc. 




FIG. 11. (Golor online). The second-order rogue wave so¬ 
lution (l30l) for (a = 2, 7 = 3, and (5 = — 1 in the vicinity of 
X = 0. (a): The motion of peaks and depressions (red and 
green lines, respectively), (b) The evolution of the absolute 
values of the peaks and depressions (red and green lines, re¬ 
spectively) . 



there are two peaks and three depressions in the inter¬ 
vals of —0.152 < X < —0.001 and 0.001 < x < 0.152, and 
there are three peaks and four depressions at \x\ < 0 . 001 . 


FIG. 13. (a) The interaction of two W-shaped solitons (|3^ : 
(b) 2D density contour plots; and (c) cross sections for two- 
solitons ([32]), for (a = (5 = 1, and 7 = —1/12. 


2. Rational soliton pair with 7 = —1/12 

For 7 = —1/12, Eq. (j30|) yields rational two-soliton 
solutions of Eq. ([2]) as 


^p2'y{x,t) = 


1 + 


G2j{x,t) + iK2j{x,t) 


D2^{x,t) 


Jxl2 


(32) 


V. CONCLUSIONS AND DISCUSSIONS 

We have explicitly obtained the first- and second- 
order RW (rogue-wave) and rational-soliton solutions of 
the integrable fifth-order NLS equation belonging to the 
NLS hierarchy. It includes the fifth-order-dispersion and 
quintic-nonlinearity terms with three arbitrary real pa¬ 
rameters [(a, 7 , and 5 in Eq. (ED]- Several familiar in- 
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FIG. 14. (a) The interaction of two solitons (|3^ : (b) 2 D 

intensity contour plots; and (c) cross sections for two-solitons 
(|32| for a = -5, 7 = - 1 / 12 , and (5=1. 


The MI (modulation instability) of the flat CW back¬ 
ground has been investigated in the framework of the 
FONLS equation, as it drives the formation of the RWs. 
The first-order RW exhibits one peak and two mutually 
symmetric depressions, while the second-order pattern 
generates an extra peak-depression pair, in a limited in¬ 
terval of the evolution variable. The location and motion 
of the peaks and depressions have been identified, in ana¬ 
lytical and numerical forms. The structure and evolution 
of the RWs were elucidated by means of the density plots. 
Moreover, for the special case 7 = —1/12, we obtain the 
corresponding one- and two-soliton solutions of Eq. 0, 
which are rational and W-shaped, being different from 
the usual bright and dark solitons. 

We expect that the predicted results may be observed 
experimentally in various physical settings described by 
the general FONLS equation. It may be interesting to 
extend the analysis by further producing more complex 
patterns, such as those corresponding to third-order RWs 
and rational solitons, which will be studied in the future. 


tegrable models, such as the classical-NLS, Hirota, and 
LPD equations are particular cases of the FONLS equa¬ 
tion, corresponding to particular values of (a, 7, and S, 
therefore, the previously studied RW solutions of those 
equations can be obtained as particular cases of our so¬ 
lutions (p!8|) and (|3Q|) . 
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APPENDIX A 


The explicit form of the system of equations {9 (|'0i(x, t)p) /dt = 0, 9 (|'0i(x,t)p) /dx = 0}: 

(t -h 6 + 30 xS) (3 — 4— 48 txa — 240 txS — 144 x^a^ — 1440 x^a S — 3600 x^S^ 

+1728 X V + 288 x^j + 12 x^) = 0, 

-64 X + 276480 t^a xS + 2488320 ta^x^S + 12441600 ta x^S^ + 663552 ta x^'f 
+110592 toLX^-i + 3317760 + 552960 tbx^^i - 256- 5308416xV 

-1769472 X V - 221184 x V “ 12288 - 1152 tot - 5760 tb - 6912 xo? - 172800 xb^ 

+1536 t^OL + 7680 + 331776 x^a^ + 207360000 x^6^ + 768 xt^ - 9216 - 1536 xq 

-69120 xa (5 + 4608 ta x^ + 230401(5 + 691200 t^S^x + 18432 xq + 110592 xq^t^ 

+49766400x^^2^2 165888 taV + 165888000x3^^3 + 207360001^3^2 

+27648 ^2^2^ + 6635520x3^3^ = 0. 


APPENDIX B 


Functions D2, G2, and 7^2, which are parts of analytical solution (|3Q|) . are given by the following expressions: 
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D 2 {x, t) = 64t^ + 384'i;xt^ + 192(5^ + bv‘^)x^t^ + 48t^ + 256i;(35^ + bv^)x^t^ — 64(120(5 + v)xt^ 
+192(^2 + bv‘^){B‘^ + v‘^)xH‘^ + 108^2 + 9 + 4(307^2 + 57600(5^ + 5280(51; 

+139^2 _ 2725 + 64)^2 + 96(85 - 115^ - 240^i; - bv‘^)xH‘^ + + v‘^fxH 

+24(240^ + llv)xt + 192(1205^^ - lB‘^v - I2{)5v‘^ - + SBv)x^t + 64(5^ 

-\-v‘^)^x^ + 16(435^ + 14405^^1; — lSB‘^v‘^ — 480(5'(;^ — 131;'^ — 165^ + 4:SBv‘^)x^, 

G 2 {x, t) = -192t^ - 1152(5^ + v‘^)xH‘^ - 76Svxt^ - 2SSt‘^ - 76Sv{3B^ + v‘^)x^t - 192(120(5 + 7v)xt 
-192(55^ + ^ 2)(52 ^2)^4 _ 240^1; + 11^^ - SB)x^ + 36, 

K2{x^ t) = —384xt^ — 7685(5^ + 3v‘^)x^t^ + 192(75 — A)xt^ — lb3bBv{B‘^ + v^)x^t — lb3bBvx‘^t^ 
-384(1205(5 - 3Bv + Av)x‘^t - 192(55^ + 2405^i; + Bv‘^ - 4B‘^ + 4v‘^)x^ 

-3845(52 ^2)2^5 24(235 - 8). 

APPENDIX C 


Functions D 2 ^{x,t)^ G 2 'y{x,t), and A 2 ^(x,t), which are parts of analytical solution (j32j), are given by the following 
expressions: 

D 2 j{x, t) = 64+ 384 xvt^ + 960 v‘^t^x^ + 48+ 1280 v^t^x^ — (7680 (5 + 64 n) t^x 
+960 vH‘^x^ — (23040 (5 n + 480 v‘^) t^x‘^ + 108+ 384 v^tx^ 

- (23040 (5 ^2 + 576 tx^ + (5760 (5 + 408 v) tx + 64 v^x^ 

- (7680 ^+ 208 v^) x^ + (256 + 230400 + 21120 (5n + 556 v‘^) x‘^ + 9, 

G 2 ^(x, t)= 36 - 192 - 768 t^vx - 1152 vH‘^x‘^ - 288 - 768 vHx‘^ - 23040 tx5 

- 1344 xvt — 192 v^x^ — 23040 x‘^d v — 1056 x‘^v‘^ ^ 

K2^{x^ t) = —192 X [4(t + nt)2 + l] . 
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